Tangling clustering of inertial particles in stably stratified turbulence 
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We have predicted theoretically and detected in laboratory experiments a new type of particle 
clustering (namely, tangling clustering of inertial particles) in a stably stratified turbulence with 
imposed mean vertical temperature gradient. In the stratified turbulence a spatial distribution of the 
mean particle number density is nonuniform due to the phenomenon of turbulent thermal diffusion, 
i.e., the inertial particles are accumulated in the vicinity of the minimum of the mean temperature of 
the surrounding fluid, and a non-zero gradient of the mean particle number density, VAT, is formed. 
It causes generation of fluctuations of the particle number density by tangling of the large-scale 
gradient, VJV, by velocity fluctuations. In addition, the mean temperature gradient, VT, produces 
the temperature fluctuations by tangling of the large-scale gradient, VT, by velocity fluctuations. 
The anisotropic temperature fluctuations contribute to the two-point correlation function of the 
divergence of the particle velocity field, i.e., they increase the rate of formation of the particle 
clusters in small scales. We have demonstrated that in the laboratory stratified turbulence this 
tangling clustering is much more effective than a pure inertial clustering (preferential concentration) 
that has been observed in isothermal turbulence. In particular, in our experiments in oscillating grid 
isothermal turbulence in air without imposed mean temperature gradient, the inertial clustering is 
very weak for solid particles with the diameter ~s 10 /im and Reynolds numbers based on turbulent 
length scale and rms velocity, Re = 250. In the experiments the correlation function for the inertial 
clustering in isothermal turbulence is much smaller than that for the tangling clustering in non- 
isothermal turbulence. The size of the tangling clusters is of the order of several Kolmogorov length 
scales. The clustering described in our study is found for inertial particles with small Stokes numbers 
and with the material density that is much larger than the fluid density. Our theoretical predictions 
are in a good agreement with the obtained experimental results. 

PACS numbers: 47.27.tb, 47.27.T-, 47.55.Hd 



I. INTRODUCTION 



Laboratory experiments [lj-|3j], numerical simulations 
0-0| and atmospheric observations [cj-fl4| revealed small- 
scale long-living inhomogeneities (clusters) in spatial dis- 
tribution of particles in different turbulent flows. The ori- 
gin of these inhomogeneities is not always clear but their 
influence on particle transport and mixing can be hardly 
overestimated. In particular, turbulence causes forma- 
tion of small-scale droplet inhomogeneities, increases the 
relative droplet velocity and affects the hydrodynamic 
droplet interactions. All these effects can enhance the 
rate of droplet collisions (see reviews [ll], HH). It is 
known that atmospheric clouds are regions with strong 
turbulence, and the turbulence effects are of a great im- 
portance for understanding of rain formation in atmo- 
spheric clouds, e.g., these effects can cause the droplet 
spectrum broade ning and acceleration of raindrop for- 
mation i, m m 0jf . 
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Different kinds of particle clustering, i.e., the preferen- 
tial concentration of inertial particles have been studied 
in a number of numerical simulations 1 1 6l— l30l j . laboratory 
experiments (see review [3l| and [32l - l38j ]) and analytical 
investigations [39l-[53| . 

Inertial clustering in isothermal and nearly isotropic 
turbulence with the Reynolds numbers Re = u$£q/v ~ 
10 3 produced by fans at the corners of the box has been 
observed in laboratory experiments (36j with hollow glass 
spheres with a mean diameter of 6 /im in air. Here Uq 
is the characteristic turbulent velocity in the maximum 
scale of turbulent motions £q and v is the kinematic vis- 
cosity. The three-dimensional radial distribution func- 
tion (RDF) as a statistical measure of clustering, has 
been determined from the particle position field using 
three-dimensional digital holographic particle imaging. 
These experiments reveled inertial particle clustering at 
the dissipation scales. The detailed comparisons between 
experiments and direct numerical simulations of iner- 
tial particle clustering in (36j have demonstrated a good 
agreement. 

The quantitative measurements of inertial clustering 
have been also performed in [33|, 38]. In particular, the 
two-dimensional RDF of particles suspended in a tur- 
bulence flow with the Reynolds numbers Re ~ 10 3 has 
been measured in [33[ by shining a laser sheet at the par- 
ticles (the glass particles with the sizes 20 /xm and 50 fim, 
the lycopodium particles with the size 25 /im), whereby 
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particle locations have been determined by a CCD cam- 
era. The experiments have been conducted in the spher- 
ical turbulence chamber with eight synthetic jet actua- 
tors which create homogeneous and isotropic turbulence 
with no mean flow at the center of the chamber. The 
experiments conducted in [33| have shown that the iner- 
tial clustering is most effective for particles with Stokes 
numbers near unity and the size of the clusters is of the 
order of ten Kolmogorov length scales. 

A one-dimensional RDF has been measured in [38| by 
sampling droplet arrivals at a fixed volume in a wind tun- 
nel, whereby the arrival statistics of water droplets with 
the mean size 22 /Ltm in turbulence with the Reynolds 
numbers Re ~ f 4 has been used to determine the one- 
dimensional RDF. In the experiments described in [38| a 
phase Doppler interferometer downstream of the active 
grid and a spray system have been used, and strong iner- 
tial clustering has been observed. Remarkably, the sim- 
ilar experimental findings have been reported in (33l . |36| 
using completely different experimental set-ups. 

The mechanism of inertial clustering of particles in 
isothermal turbulence is as follows [54[ . The particles in- 
side the turbulent eddies are carried out to the boundary 
regions between the eddies by the inertial forces. This 
mechanism of the inertial clustering acts in all scales of 
turbulence, and is more pronounced in small scales. 

The goal of this study is to investigate experimen- 
tally and theoretically a new type of particle clustering, 
namely tangling clustering of inertial particles in stably 
stratified turbulence with imposed mean vertical temper- 
ature gradient. In experimental study of particle cluster- 
ing we use Particle Image Velocimetry to determine the 
turbulent velocity field, an image processing technique 
to determine the spatial distribution of particles, and a 
specially designed temperature probe with twelve sensi- 
tive thermocouples for measurements of the temperature 
field. The clustering described in our study is found for 
inertial particles with small Stokes numbers and with ma- 
terial density that is much larger than the fluid density. 

In the stratified turbulence a spatial distribution of 
the mean particle number density is nonuniform due to 
phenomenon of turbulent thermal diffusion [55l |56| . In 
particular, the inertial particles are accumulated in the 
vicinity of the minimum of the mean temperature of the 
surrounding fluid, which causes formation a non-zero gra- 
dient of the mean particle number density, VAT. The 
phenomenon of turbulent thermal diffusion has been pre- 
dicted theoretically in [55|, detected in the laboratory 
experiments in stably and unstably stratified turbulent 
flows in [57H60| and observed in atmospheric turbulence 
in [6l|. 

Fluctuations of the particle number density can be gen- 
erated by tangling of the gradient, VA, of the mean 
particle number density by velocity fluctuations • On 
the other hand, the imposed mean temperature gradient, 
VT, results in generation of the anisotropic temperature 
fluctuations by tangling of this large-scale gradient, VT, 
by velocity fluctuations. These temperature fluctuations 



may contribute to the two-point correlation function of 
the divergence of the particle velocity field. The latter 
enhances the rate of formation of particle clusters in small 
scales by the tangling mechanism. 

The tangling mechanism is universal and independent 
of the way of generation of turbulence for large Reynolds 
numbers. For instance, tangling of the gradient of the 
large-scale velocity shear produces anisotropic velocity 
fluctuations [63|, [64[ which are responsible for different 
phenomena: formation of large-scale coherent structures 
in a turbulent convection |65|, excitation of the large- 
scale inertial waves in a rotating inhomogeneous tur- 
bulence [gfjj], generation of large-scale vorticity [13] and 
large-scale magnetic field [g8[ in a sheared turbulence. 

The paper is organized as follows. Section II describes 
the experimental set-up for a laboratory study of the tan- 
gling clustering of inertial particles in stably stratified 
turbulence. The data processing and experimental re- 
sults are presented in Section III. The theoretical analysis 
and comparison with experimental results are performed 
in Section IV. Finally, conclusions are drawn in Section 
V. 



II. EXPERIMENTAL SET-UP 

The experiments were carried out in a turbulence gen- 
erated by oscillating grids in air. The test section of the 
oscillating grids turbulence generator was constructed as 
a rectangular chamber with dimensions 29 x 58 x 29 cm 3 
(see Fig. [Tj). Pairs of vertically oriented grids with bars 
arranged in a square array (with a mesh size 5 cm) are 
attached to the right and left horizontal rods driven by 
speed-controlled motors. The grids are positioned at a 
distance of two grid meshes from the chamber walls par- 
allel to them. Both grids are operated at the same am- 
plitude of 61 mm, at a random phase and at the same 
frequency varied in the range from 2.2 Hz to 16.5 Hz. 
Here we use the following system of coordinates: Z is 
the vertical axis, the Y"-axis is perpendicular to the grids 
and the XZ-plane is parallel to the grids. 

A mean temperature gradient in the turbulent flow was 
formed with two aluminium heat exchangers attached to 
the bottom and top walls of the chamber. We performed 
experiments in stably stratified fluid flow (the cold bot- 
tom and hot top walls of the chamber). In order to im- 
prove heat transfer in the boundary layers at the walls we 
used two heat exchangers with rectangular fins (3x3x15 
mm , see Fig. 1) which affords a mean temperature gra- 
dient 118 K/m at a mean temperature of about 308 K. 
All experiments were conducted at the same temperature 
difference between the top and bottom walls AT = 50 K 
(e.g., the bottom wall temperature was 283 K and the 
top wall temperature was 333 K). 

The temperature field was measured with a tempera- 
ture probe equipped with twelve E-thermocouples (with 
the diameter of 0.13 mm and the sensitivity of 65 /iV/K) 
attached to a vertical rod with a diameter 4 mm. The 
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FIG. 1: Experimental set-up: (1) -laser light sheet; (2) - heat 
exchangers; (3) - grid driver; (4) - digital CCD camera. 



spacing between thermocouples along the rod was 22 mm. 
Each thermocouple was inserted into a 1 mm diameter 
and 45 mm long case. A tip of a thermocouple protruded 
at the length of 15 mm out of the case. The mean tem- 
perature was measured for 5 rod positions with 50 mm 
intervals in the horizontal direction, i.e., at 60 locations 
in a flow. The exact position of each thermocouple was 
measured using images captured with the optical sys- 
tem employed in PIV measurements. A sequence of 1024 
temperature readings (each reading was averaged over 15 
instantaneous measurements) for every thermocouple at 
every rod position was recorded and processed using the 
developed software based on Lab VIEW 7.0. 

The turbulent velocity field was measured using a dig- 
ital Particle Image Velocimetry (PIV) system with LaVi- 
sion Flow Master III (see, e.g., |69l - [7ll |). A double-pulsed 
Nd-YAG laser (Continuum Surelite 2 x 170 mJ) is used 
for light sheet formation. Light sheet optics comprise 
spherical and cylindrical Galilei telescopes with tune- 
able divergence and adjustable focus length. We em- 
ployed a progressive-scan 12 Bit digital CCD camera 
(pixels with a size 6.7 /zm x 6.7 /an each) with dual frame 
technique for cross-correlation processing of captured im- 
ages. The tracer used for PIV measurements was incense 
smoke with sub-micron particles (with the material den- 
sity ptr ~ 1 g/cm 3 ), which was produced by high tem- 
perature sublimation of solid incense particles. Analysis 
of smoke particles using a microscope (Nikon, Epiphot 
with an amplification 560) and PM-300 portable laser 
particulate analyzer showed that these particles have an 
approximately spherical shape with a mean diameter of 
0.7/im. 

We determined the mean and the r.m.s. velocities, two- 
point correlation functions and an integral scale of tur- 
bulence from the measured velocity fields. Series of 130 
pairs of images, acquired with a frequency of 4 Hz, were 
stored for calculating velocity maps and for ensemble and 



spatial averaging of turbulence characteristics. The cen- 
ter of the measurement region coincides with the center 
of the chamber. We measured velocity for flow areas from 
50 x 50 mm 2 up to 212 x 212 mm 2 with a spatial resolution 
of 1024 x 1024 pixels. This corresponds to a spatial res- 
olution from 48 fim / pixel up to 207 /im / pixel. These 
measurement regions were analyzed with interrogation 
windows of 32 x 32 or 16 x 16 pixels, respectively. 

In every interrogation window a velocity vector was 
determined from which velocity maps comprising 32 x 32 
or 64 x 64 vectors were constructed. The mean and r.m.s. 
velocities for every point of a velocity map (1024 points) 
were calculated by averaging over 130 independent maps, 
and then they were averaged over 1024 points. The two- 
point correlation functions of the velocity field were cal- 
culated for every point of the central part of the velocity 
map (with 16 x 16 vectors) by averaging over 130 in- 
dependent velocity maps, and then they were averaged 
over 256 points. An integral scale £q of turbulence was 
determined from the two-point correlation functions of 
the velocity field. 

Particle spatial distribution was determined using dig- 
ital Particle Image Velocimetry (PIV) system. In partic- 
ular, the effect of Mie light scattering by particles was 
used to determine the particle spatial distribution in the 
flow (see, e.g., [72}). In the experiments we probed the 
central 212 x 212 mm 2 region in the chamber. The mean 
intensity of scattered light was determined in 32 x 16 
interrogation windows with the size 32 x 64 pixels. The 
vertical distribution of the intensity of the scattered light 
was determined in 16 vertical strips composed of 32 in- 
terrogation windows. The light radiation energy flux 
scattered by small particles is E s oc E ^(Trd p /X;ao;N), 
where Eq oc 7rei 2 /4 is the energy flux incident at the par- 
ticle, dp is the particle diameter, A is the wavelength, clq 
is the index of refraction and *Sf is the scattering func- 
tion. For wavelengths A which are larger than the par- 
ticle perimeter (A > 7rc? p ), the function ^ is given by 
Rayleigh's law, ^ oc dp. If the wavelength is small, the 
function ^ tends to be independent of d p and A. In the 
general case the function ^ is given by Mie's equations 
(sec, e.g., 73], Chapter 4). The scattered light energy 
flux incident on the CCD camera probe (producing pro- 
portional charge in every CCD pixel) is proportional to 
the particle number density n, i.e., E s cx Eq n (irdp/4). 

In order to characterize the spatial distribution of par- 
ticle number density n cx E T /E in the non-isothermal 
flow, the distribution of the scattered light intensity E 
measured in the isothermal case was used for the nor- 
malization of the scattered light intensity E T obtained 
in a non-isothermal flow under the same conditions. The 
scattered light intensities E T and E in each experiment 
were normalized by corresponding scattered light intensi- 
ties averaged over the vertical coordinate. Mie scattering 
is not affected by temperature change because it depends 
on the electric permittivity of particles, the particle size 
and the laser light wave length. The temperature effect 
on these characteristics is negligibly small. Similar exper- 
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imental set-up and data processing procedure were used 
in experimental study of different aspects of turbulent 
convection [zH LZ5| and in [57T - l60l ] for investigating the 
phenomenon of turbulent thermal diffusion [5a . l56j . 

For experimental study of particle clustering we used 
hollow borosilicate glass particles having an approxi- 
mately spherical shape, a mean diameter of 10 /im and 
the material density p p w 1.4 g/cm 3 . These particles 
have been injected in the chamber using an air jet in or- 
der to improve particle mixing and prevent from particle 
agglomeration. 



III. DATA PROCESSING AND 
EXPERIMENTAL RESULTS 



All experiments for study of the particle clustering 
have been performed at the frequency of grid oscillations 
/ = 10.4 Hz. The turbulent flow parameters in the os- 
cillating grids turbulence generator at / = 10.4 Hz are 
as follows: the r.m.s. velocity is uq = 

VTu 5 ! = 12 cm/s, 

the integral (maximum) scale of turbulence is £q = 3.2 
cm, the Reynolds numbers Re = Uq£q/v = 250, the Kol- 
mogorov length scale is £ n = £o/Re 3 ^ 4 — 510 /im and the 
Kolmogorov time scale r n = To/Re 1 / 2 = 1.7 x 10~ 2 s, 
where r = (q/ u o- The Stokes time for the particles with 
the diameter d p = 10 /im is t s = 10~ 3 s, the Stokes num- 
ber St = t s /t v = 5.9 x 10 -2 , the coefficient of molecular 
diffusion D m = 1.4 x 10~ 8 cm 2 /s and the Peclet number 
Pc= u £o/D m = 3 x 10 9 . 

The velocity measurements have shown that there is a 
slight difference in the velocity components (about 9%) 
caused by an anisotropy of forcing of the grid oscillating 
turbulence. The inhomogeneity of turbulence in the core 
of fluid flow is weak. We have found a weak mean flow 
in the form of two large toroidal structures parallel and 
adjacent to the grids. The interaction of these structures 
results in a symmetric mean flow that is sensitive to the 
parameters of grids adjustment. We studied the parame- 
ters that affect the mean flow, e.g., the grids distance to 
the walls of the chamber, the angles of the grids planes 
with the axes of their oscillation. This study allowed us to 
expand the central region with homogeneous turbulence 
by inserting partitions behind the grids. The measured 
r.m.s. velocity was 5 times higher than the characteristic 
mean velocity in the core of the flow. 

For analysis of particle clustering we use the 
radial distribution function (RDF), G(R) = 
(n(t,x)n(t,y))/N(t,x)N(t,y), that is the condi- 
tional probability density of finding a second particle 
at a give n separation distance from a test particle (see, 
e.g., [76[), where n(t, x) is the instantaneous number 
density of particles, N(t,x) = (n(t,x)) is the mean 
number density of particles, the angular brackets denote 
ensemble averaging and R = y x. The RDF can be 
estimated from a field of M particles by binning the 
particle pairs according to their separation distance, so 



that the function G(R) is estimated as follows: 



G(R) 



N^/V 

r(p) 



(1) 



(see, e.g., .77]), where N^ v is the number of particle 
pairs separated by a distance R ± AR/2, AV is the 
volume of the spherical shell located between i?± AR/2, 
= M(M - l)/2 is the total number of pairs and V 
is the total volume of the probed region. 

In our experiments we employ the PIV system in order 
to determine the particle spatial distribution. Since we 
use the two-dimensional images, Eq. (fT]) is modified as 
follows: 



G(R) 



N^l/AS 
N { s p) /S 



(2) 



where AS" = tt[{R+ AR/2) 2 - (R- AR/2) 2 } is the area of 
the annular domain located between i?± AR/2, S is the 
area of the part of the image with the radius R max that 
was used in data processing in order to exclude the edge 
effects. In our experiments with the image sizes 20 x 20 
cm 2 and 5x5 cm 2 , the maximum radius i? ma x = 0.8 
cm and i? m ax = 0.2 cm, respectively. The total number 
of particles in the image 5x5 cm 2 is of the order of 
M ~2x 10 4 . Using the thickness of the laser light sheet 
d = 0.2 cm we can estimate the effective 3D particle mean 
number density in the experiments as N ~ 4 x 10 3 cm~ 3 . 

The measured radial distribution function has been 
used for determining the two-point correlation function of 
the particle number density, $(£, R) = (0(t, x)0(t, x + 
R)), where 0(t, x) = n(t,x) — 7V(t,x) is the deviation 
of the instantaneous number density of particles n(t, x) 
from the mean number density of particles iV(t, x). The 
two-point correlation function of the particle number 
density is given by 



$(t,R) = iV 2 [G(i,R)-l] 



(3) 



(see, e.g., [76|). Note that in our experiments the correla- 
tion function has been normalized by squared mean num- 
ber density of particles in every image. We perform the 
double averaging over all particles in the image and then 
over ensemble of 50 images, which allow us to increase the 
accuracy of determining the two-point correlation func- 
tion of the particle number density. Therefore, particle 
clustering in our study is understood in the statistical 
sense by applying an ensemble averaging over many im- 
ages with the instantaneous particle distributions. 

Since a typical size of a particle cluster is of the order 
of several Kolmogorov length-scales of turbulence (see, 
e.g., [H,|3l,|43]), we have to use a sub- pixel resolution in 
the data analysis. In our experiments one pixel is of the 
order of 1/3 of the Kolmogorov length scale of turbulence. 
On the other hand, the size of the analyzed region in the 
image cannot be reduced strongly, because a number of 
particles in the analyzed region of the image should be 
large in order to provide a good statistics. 



FIG. 2: Normalized two-point second-order correlation func- 
tion <&(R) determined in our experiments (filled squares) and 
the least-square fit for the experimental results (dashed line). 
Normalized two-point second-order correlation function $(i?) 
determined from our model for a T = 1, a c — 1 (solid line). 
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FIG. 3: Normalized two-point second-order correlation func- 
tion <&(R) determined in our experiments: (i) for isothermal 
turbulence (unfilled circles); (ii) for non- isothermal turbulence 
(filled squares as in Fig. [2]). 



In order to attain the sub-pixel resolution we employ 
the following method: (i) we determine the response 
function for the CCD camera used in the PIV system 
by analyzing the light intensity distribution in the im- 
age for single particles located at the center of the pixel 
in the form of the Gaussian distribution with er» = 0.81 
pixels; (ii) segmentation of the image using a threshold 
technique; (iii) identification of particle locations in the 
segments by least-square fitting of the recorded light in- 
tensity distribution and the light intensity distribution 
caused by superposition of the Gaussian distributions at 
the particle locations. This procedure has been tested 
using artificial computer-generated images. The error in 
the determining the particle coordinates is of the order 
of 3%. 

This approach allows us to determine in the exper- 
iments the two-point second-order correlation function 
$(i?) of the particle number density. For instance, in 
Fig. [2] we plot the normalized two-point second-order cor- 
relation function $(-/?) determined in our experiments, 
which arc performed in the air flow with imposed mean 
temperature gradient, i.e., for the stably stratified tur- 
bulence with the temperature difference, AT = 50 K, 
between the top and bottom walls of the chamber. 

On the other hand, we have found that in the exper- 
iments in isothermal turbulence without imposed mean 
temperature gradient, the inertial particle clustering is 
very weak, i.e., the correlation function $(-R) for the iner- 
tial clustering is much smaller than that for the tangling 
clustering. Indeed, in Fig. [3] we compare the correlation 
function $(i?) for the isothermal turbulence (unfilled cir- 
cles) with that for the turbulence with imposed mean 
temperature gradient (filled squares). For instance, for 
the isothermal turbulence <&i n (R — 0A£ n ) = 0.3, while 
for the turbulence with imposed mean temperature gra- 
dient $tan(-R = 0.44,) = 3.3 (the tangling clustering). 
The minimum distance Rq at which the correlation func- 



tion $(i?) approaches zero is Rq = £ v for the isothermal 
turbulence and Rq — 7.6£ v for the tangling clustering. 

In the next section we perform theoretical study of 
this new kind of particle tangling clustering and compare 
theoretical predictions with the experimental results. 



IV. THEORETICAL ANALYSIS AND 
COMPARISON WITH EXPERIMENTAL 
RESULTS 

Let us consider the two-point second-order correlation 
function of the particle number density fluctuations gen- 
erated by tangling of the gradient of the mean particle 
number density by the turbulent velocity field. This gra- 
dient is formed due to phenomenon of turbulent thermal 
diffusion in stably stratified turbulence. 



A. General consideration 

The equation for the number density n(t, x) of particles 
advected by a fluid velocity field reads: 

d ti 

+ V-(nv) =D m An, (4) 

at 

where D m is the coefficient of molecular (Brownian) dif- 
fusion and v(i,x) is the particle velocity field. Equa- 
tion (jj) implies conservation of the total number of par- 
ticles in a closed volume. The equation for fluctuations of 
the particle number density, 0(i,x) = n(t, x) — N(t, x), 
reads: 

^ + V-(6v-(ev))=£) m Ae 

-(v-V)iV-iVV-v, (5) 
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(see, e.g., (zl, [79j]). Using Eq. ([5]) we derive equation 
for the evolution of the two-point second-order correla- 
tion function of the particle number density, R) = 
(e(t,x)0(f,x + R)): 



~dt 



[B(R) + 2U^(R) ■ V + AjfRjV.Vj] $(i,R) 
+/(R) , (6) 
(see 0), where U^(R) = (1/2) [U(R) - U(-R)], 

Dij = 2D m * tf +Dj(0)-Dj(R), (7) 

/>GO 

Dj (R) = 2 / ( Wi [0, £(t, x|0)] Wj [r, €(t, x + R|r)] ) dr, 
Jo 

(8) 

/•OO 

R(R) = 2 / (6[0,€(t,x|0)]6[r,€(t,x + R|r)])dr, 
Jo 

(9) 

POO 

Ui(R) = -2 / (« i [0,€(t,x|0)]6[T,€(t,x + R|T)])dT. 
Jo 

(10) 

T 

6 = divv, Z).-(R) is the scale-dependent turbulent dif- 
fusion tensor, Sij is the Kronecker tensor, /(R) is the 
source of particle number density fluctuations and (...) 
denotes averaging over the statistics of turbulent velocity 
field and the Wiener process w(£). The Wiener trajec- 
tory £(i,x|s) (which is often called the Wiener path) in 
the expressions for the turbulent diffusion tensor (R) 
and other transport coefficients is defined as follows: 

£(t,x|«) = x- J v[t,€(*,x|t)] dr - ^2D~w(t- s), 

(11) 

where w(i) is the Wiener random process which describes 
the Brownian motion (molecular diffusion). The Wiener 
random process w(i) is defined by the following proper- 
ties: (w(£)) w = , (wi(t + T)w 3 (t)) w = r<%, and (. . .) w 
denotes the mathematical expectation over the statistics 
of the Wiener process. The velocity Vi[r, x|r)] de- 
scribes the Eulerian velocity calculated at the Wiener 
trajectory. 

The source function /(R) in Eq. © is related to the 
last two terms — (v-V)iV — 7VV- v in the right hand side 
of Eq. (|5j|. In particular, when ViV ^ the nonzero 
source I results in generation of fluctuations of the par- 
ticle number density caused by tangling of the gradient 
of the mean particle number density by the turbulent 
velocity field. The source function /(R) is given by: 



J(R) = B(R)N 2 +U (S) (R) ■ ViV 2 
+§I>J(R) (V;iV) (VjN), 



(12) 



where U^(R) = (1/2) [U(R) + U(-R)] and we 
have taken into account that V, (x) vj y) iV(i, x)N(t, y) = 
(3/4) (VjiV) (V,W) and V, = vj R) . 



The meaning of the turbulent transport coefficients 
5(R) and U(R) is as follows. The function B(R) is 
determined by the compressibility of the particle veloc- 
ity field. The vector U(R) determines a scale-dependent 
drift velocity which describes transport of fluctuations 
of particle number density from smaller scales to larger 
scales, i.e., in the regions with larger turbulent diffusion. 
The scale-dependent tensor of turbulent diffusion (R) 
in very small scales is equal to the tensor of the molecular 
(Brownian) diffusion, while in the vicinity of the maxi- 
mum scale of turbulent motions this tensor coincides with 
the regular tensor of turbulent diffusion. 

If the turbulent velocity field is not delta-correlated 
in time (e.g., the correlation time is small yet finite), 
the tensor of turbulent diffusion, D,^ (R) , is compress- 
ible, i.e., (d/dRi)D^{R) ^ (for details see The 
parameter a T that characterizes the degree of compress- 
ibility of the tensor of turbulent diffusion is defined as 
follows: 



V,V,/> ,.R 



VjVjD mn (R)e imp e 



J np 



(13) 



where is the fully antisymmetric Levi-Civita unit ten- 
sor, £ = £ x with \t — s\ t c . Here t c is the charac- 
teristic turbulent time, 



5«,[r,£(i,x|r)] d£ v 



dxi 



dr 



and 



dxi 



r* du m [r',HtMr')] <Kn dT , 
J T d^n dxi 



(14) 



(15) 



When the turbulent velocity field is a delta-correlated in 
time random process, V-£ = — (V-u) (t — s) and Vx£ = 
— (Vxu) (t — s). Therefore, in this case the parameter 
a T is given by 



((V-u) 2 ) 

= (Tii 



((Vxu) 2 ) 



(16) 



where u is the fluid velocity field. For a small yet finite 
correlation time (i.e., for a small Strouhal numbers, Sr = 



T c y/ (u 2 ) /£ <C 1), the parameter a T can be estimated as 



OV = (T, u 



For derivation of Eq. ([TT]) we used Eq. (C12) given in 
42]. For inertial particles with small Stokes time the 
corrections ~ 0(St 2 ) in Eqs. (|T4| - (fT6|) can be neglected 
(see subsection IV-B), where St = t s /t v is the Stokes 
number, t s is the Stokes time of inertial particles and T n 
is the Kolmogorov time scale. 

Equation © with /(R) = has been derived in [8(| for 
a delta-correlated in time random incompressible (b = 0) 
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velocity field. For a turbulent compressible (6 ^ 0) veloc- 
ity field with a finite correlation time Eq. ([5]) has been de- 
rived in [13] by means of stochastic calculus, i.e., Wiener 
path integral representation of the solution of the Cauchy 
problem for Eq. using Feynman-Kac formula and 
Cameron-Martin-Girsanov theorem. The comprehensive 
description of this approach can be found in [42|, [62|, [8l[ . 



B. Gradient of the mean particle number density 

A nonzero gradient of the mean particle number den- 
sity in stratified turbulence with external mean temper- 
ature gradient is caused by the phenomenon of turbulent 
thermal diffusion [55l - [6lj | . This phenomenon in turbu- 
lent stratified flows results in the non-diffusive flux of 
particles in the direction of the heat flux. Particles are 
accumulated in the vicinity of the minimum of the mean 
temperature of the surrounding fluid. This causes for- 
mation of large-scale inhomogeneities in spatial distribu- 
tion of particles. Turbulent thermal diffusion has been 
detected in two experimental set-ups: oscillating-grids 
turbulence generator (57l - [59| and multi-fan turbulence 
generator [60] . The experiments have been performed 
for stably and unstably stratified fluid. In these exper- 
iments, even with strongly inhomogeneous temperature 
fields, particles in turbulent fluid accumulate in the re- 
gions of temperature minima, in a good agreement with 
the theory of turbulent thermal diffusion [55|, [56| . 

Let us discuss the physics of the phenomenon of tur- 
bulent thermal diffusion. The velocity of particles, v, 
depends on the velocity of the surrounding fluid, u, and 
it can be determined from the equation of motion for a 
particle. When p p 3> p, this equation represents a bal- 
ance of particle inertia with the fluid drag force produced 
by the motion of the particle relative to the surrounding 
fluid, dv/dt = (u — v)/t s , where t s is the particle Stokes 
time, p is the fluid density and p p is the material density 
of a particle. Solution of the equation of motion for small 
Stokes numbers, St <C 1, reads: 



= u-St 



- + (u.V)u 



0(St 2 ) 



(18) 



(see, e.g., 54]). This solution is written in dimension- 
less form, where the time is measured in the units of 
Kolmogorov time scales. The second term in Eq. ([18]) 
describes the difference between the local fluid velocity 
and particle velocity arising due to the small but finite in- 
ertia of the particle. In this study we consider low Mach 
numbers turbulent flow with V u = — p^ 1 (u-V)p ^ 0. 



Equation ([T5|) for the velocity of particles and Navier- 
Stokes equation for the fluid for large Reynolds numbers 
yield the equation for V- v written in dimensional form: 



V-v = 



V- u r s V- 
1 



(u-V)p+-V z p + 0{ri) 



(19) 



where p is the fluid pressure. 

The physical mechanism of the phenomenon of tur- 
bulent thermal diffusion for inertial particles can be ex- 
plained as follows. Due to inertia, particles inside the tur- 
bulent eddies drift out to the boundary regions between 
the eddies (the regions with the decreased velocity of the 
turbulent fluid flow). Neglecting non-stationarity and 
molecular viscosity, the estimate based on the Bernoulli's 
law implies that these are the regions with the increased 
pressure of the surrounding fluid. Consequently, particles 
are accumulated in the regions with the maximum pres- 
sure of the turbulent fluid. Indeed, due to the inertia ef- 
fect V- v cx (r s /p) V 2 p 7^ even for incompressible fluid 
flow [see Eq. ([T9l ]. On the other hand, for large Peclet 
numbers, when we can neglect the molecular diffusion of 
particles in Eq. V- v oc —dn/dt. This implies that in 
regions with maximum pressure of turbulent fluid (i.e., 
where V 2 jj < 0) there is accumulation of inertial parti- 
cles (i.e., dn/dt cx — (r s /p)V 2 p > 0). Similarly, there 
is an outflow of inertial particles from regions with the 
minimum pressure of fluid. 

In case of homogeneous and isotropic turbulence with- 
out external large-scale gradients of temperature, a drift 
from regions with increased (decreased) concentration of 
particles by a turbulent flow of fluid is equiprobable in 
all directions. Therefore pressure (temperature) of the 
fluid is not correlated with the turbulent velocity field 
and there exists only turbulent diffusion of particles. 

Situation drastically changes in a turbulent fluid with 
a mean temperature gradient. In this case, the heat flux 
(u9) is not zero, i.e., fluctuations of fluid temperature, 
8, and velocity are correlated. We consider low-Mach- 
number flows (M = u/c s <C 1, c s is the sound speed) and 
study mean- field effects. For low-Mach- number isother- 
mal flows, the mean fluid mass flux (u p') is very small 
(~ 0(M 2 )) (see, e.g., [8a|). i.e., the fluctuations of 
the fluid density p' and velocity u are weakly corre- 
lated. Moreover, in stratified turbulent flows with im- 
posed mean temperature gradient the total mass current 
in the cell reference frame vanishes. 

On the other hand, fluctuations of pressure must be 
correlated with the fluctuations of velocity due to a 
non-zero turbulent heat flux, (u(9) ^ 0. Indeed, us- 
ing the equation of state for an ideal gas we find that 
p/P = p'/p + e/T, and (up)/P = (u0)/T, where P, 
T and p are the mean fluid pressure, temperature and 
density, respectively. Therefore, the fluctuations of tem- 
perature and pressure are correlated and the pressure 
fluctuations cause fluctuations of the number density of 
particles. The correlation between p and 9, necessary for 
P(9u) = T(pu), arises from the buoyancy component of 
p and from the effect of non-uniform mass density in the 
Navier-Stokes equation. 

Increase (decrease) of the pressure of surrounding fluid 
is accompanied by accumulation (outflow) of the par- 
ticles, respectively. The direction of the mean flux of 
particles coincides with the direction of the heat flux of 
temperature - towards the minimum of the mean tern- 
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perature. Therefore, the particles are accumulated in 
this region (for more details, see (55j). 

Equation for the evolution of the mean number density 
N of particles reads: 



dN 
~dt 



V- [7VV + F (ti) ] = 



(20) 



where V is the mean particle velocity, = (v 0) is the 
turbulent flux of particles that includes contributions of 
turbulent thermal diffusion and turbulent diffusion (see 

MM)- 



p (n) = yeff N _ D VN 



D T ss £ a u is the coefficient of turbulent diffusion, u 
is the characteristic turbulent velocity in the maximum 
scale £q of turbulent motions, V is the effective veloc- 
ity caused by turbulent thermal diffusion is given by the 
following equation: 



V 



elf 



(V-v)) 



(22) 



Equation (|22|) for the effective velocity has been derived 
using different rigorous methods in [H, [H, HH EH H3] ■ 
Note that even a simple dimensional analysis yields the 
estimate for the effective velocity V cff that coincides with 
Eq. (HU). Indeed, the magnitude of dS/dt + V Q 
D m V 2 in Eq. © can be estimated as O/r. Therefore, 
the turbulent component of particle number density 
is of the order of w -tV-(Nv) = -r[N(V-v) + 
(v-V)TV]. Now let us determine the turbulent flux of 



particles F q 



(n) 



(Vi 0): 



F. 



(n) 



-Nr(v t (V-v)) -riviV^WjN 



(23) 



where the first term in the right hand side of Eq. (|23j) de- 
termines the turbulent flux of particles caused by turbu- 
lent thermal diffusion: -Nr(vi (V-v)) = V° s N, while 
the second term in the right hand side of Eq. (1231) deter- 
mines the turbulent flux of particles caused by turbulent 
diffusion: r (viVj)VjN = D T ViN. In the latter estimate 
we neglected the anisotropy of turbulence for simplicity. 
More detailed analysis shows that the effective velocity 
caused by turbulent thermal diffusion is given by: 



-a D n 



VT 



(24) 



where V cff = U(R = 0) (see subsection IV-C). The tur- 
bulent thermal diffusion ratio, a, is given by the following 
equation: 



a = I 



-yWgLp ln(Re) 
wo 4 



(25) 



(see [55j, |5 
heats, W 



), where 7 = c p /c v is the ratio of specific 
g — t s g is the terminal fall velocity of particles, 
m p /6irp va p is the Stokes time for small spherical 

g is the acceler- 



particles of the radius a p and mass m p , 



(21) ation of gravity, L p = \V Z P/P\ and Re = uq£q/v is the 



Reynolds number. For gases and non-inertial particles 
a = 1. For derivation of Eqs. (|2"4"|) and (|23|) we took into 
account the equation of state, neglected the mass flux 
of fluid, (up'), for the low-Mach-number flows, and used 
the identity r s = pW g Lp/P, where |V 2 P| = pg. The 
steady-state solution of Eq. ([20| at V = reads: 



ViV 



VT 



(26) 



26l) in order to determine the source 



We will use Eq. 
function I(R). 



C. Functions S(R) and U (S) (R) 

Let us consider the case a 2 ^> 1. The magnitude of a 2 
in the experiment was of the order of 10 - 15. Note that 
it is not easy to determine the exact value of a 2 in the 
experiments because of a size distribution of particles. 
In particular, the particle size in the experiments varies 
from 3 to 40 fim, with the mean value 10 pui. On the 
other hand, in the analysis of the particle clustering we 
used a threshold in the light intensity. This implies that 
we did not take into account very small particles in the 
data analysis of the particle clustering. 

In this case V-v w (r s / p)V 2 p, and the functions 
B(R) determined by Eq. (j9)), is given by the following 
equation: 



B(R) 



■It. 



■It. 



2 r p2 



, N r[V 2 p(x)] v 2 P (y)> - , 

p p z 
p2 

p 



P 2 



(t[V 2 6( X )} V 2 8(y)) + — (r[VV(x)] V 2 p\y)) 



T 2 v n w " ' p 2 



+ ^ ((^[V 2 p'(x)] V 2 8(y)) + (t[V 2 9( X )] V 2 p'(y)) 



(27) 



where V 2 p(x) = [V*- 3 ^] 2 p(x). Hereafter we omit the argument t in the correlation function. In derivation of 
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this equation we used the relationship 



P_ 
P 



EL 
p 



T 



O{p'0), 



(28) 



that follows from the equation of state for ideal gas. We 
also take into account that characteristic spatial scales for 
fluctuations of fluid pressure, temperature and density 
are much less than those for the mean fields. 

In turbulence with imposed turbulent heat flux 
(e.g., the imposed mean temperature gradient in the 
oscillating grid turbulence), the correlation function 
([V 2 6>(x)] V 2 #(y)) is much larger than the correla- 
tion functions of density-density fluctuations or density- 
temperature fluctuations, i.e., 



i |([V 2 0(x)] V 2 9{y))\ » 1 |([VV(x)] VV(y))| 



- |([V 2 0(x)] V 2 9(y))\ » 1 |([V 2 p'(x)] V 2 9(y))\ 



i |([V 2 0(x)] V 2 9{y))\ » I |([V 2 tf(x)] VV(y))| 



(29) 
(30) 
(31) 



Indeed, the correlation function ([V 2 #(x)] V 2 #(y)) is 
caused by the turbulent heat flux, i.e., (0(x) 0(y)) oc 
-T (it*(x)0(y)) (ViT) [see Eqs. (03} and (05])]. where r 
is the characteristic turbulent time. On the other hand, 
the correlation functions of density-density fluctuations 
or density-temperature fluctuations are nearly indepen- 
dent of the turbulent heat flux, and they are proportional 
to the mass flux (u(x) p'(y)), which is very small. 

In particular, the temperature fluctuations can be es- 
timated as 9 oc —To Ui VjT. Then the temperature- 
density fluctuations are estimated as (6*(x)p'(y)) oc 
—r («i(x) p'{y)) (ViT). The density fluctuations are de- 
termined by the continuity equation: 



= -V-(pu) + 0(p'u) 



(32) 



and the correlation function of density-density fluctua- 
tions (p'(x) p'(y)) is determined by the following equa- 
tion 

!<p'(x)p'(y))=-p[V^ (p'(x)u i (y)) 
+vf> (p'(y) u . i( x))]-^[<p'(x) Ul (y)) 



V(yWx))] , 



(33) 



which follows from Eq. (|32|) . Since (p'(x)Ui(y)) is 
very small, and is nearly independent of the turbu- 
lent heat flux, the correlation functions of the density- 
density fluctuations or density-temperature fluctuations 
are much smaller than the correlation functions of the 
temperature-temperature fluctuations. 



In k space the correlation function 
(t[V 2 6<(x)] [V 2 6»(y)]) reads: 



(r[V 2 0(x)] [V 2 9(y)}) = J r(k)k i (9(k)9(-k)) 

x exp (ik-R) dk. (34) 

Similarly, the function IT S )(R) determined by Eq. (jlOp . 
is given by the following equation: 



2r, 



P 

2t s P 
' PT 



ru(x) [V 2 p(y)}) 



(ru(x) [V 2 9(y)]) 



(35) 



In k space the correlation function (ru(x) [V 2 #(y)]) 
reads: 

(ru(x)[V 2 0(y)]) = - Jr(k)k 2 {u(k)9(-k)) 

x exp (ik-R) dk . (36) 

In order to determine the correlation functions 
(0{k) 9(—k)} and (u(k) 9(— k)) we use the following evo- 
lutionary equation for the temperature field Tt t(i, r) in 
a turbulent flow: 



dT to 

at 



(u- V) T tot + ( 7 - 1) (V-u) T tot - D V 2 T tot , 



(37) 



where D is the coefficient of molecular temperature dif- 
fusion, 7 = c p /c v is the ratio of specific heats, u is the 
fluid velocity field that satisfies to continuity equation in 
anelastic approximation for a low-Mach-number flow: 



V-(pu) = 



(38) 



Combining Eqs. (|3"T|) and (j38"j) we obtain the following 
equation: 



at 



(u-V)T tot = DV 2 T u 



(39) 



where u = 7U. Averaging Eq. (13T)|) over an ensemble of 
turbulent velocity field we obtain the equation for the 
evolution of the mean temperature field T(t,r): 



BT 

-77 — hV-F = D V 2 T 

at 



(40) 



where F = (u0) is the heat flux, and for simplicity we 
consider the fluid flow with a zero mean velocity. Note 
that for a low-Mach-number flow without imposed exter- 
nal pressure gradient Vp/p s» — VT/T. 

Subtracting Eq. (|4"0")) from Eq. (131)1) yields equation for 
the temperature fluctuations: 



^ + Q-DV 2 9 = I 
at 



(41) 
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where I = — (u-V)T is the source term and Q = V- [u0— 
(u(9)] + #(u-V)yo/p is the nonlinear term. In Eq. (|4"0")) 
we neglected the term (9u)-'Vp/p which is quadratic in 
large-scale spatial derivatives. 

Now we derive formulae for the function Fi(t, k) = 
(ui(t, k) 9(t, — k)) and the temperature fluctuation func- 
tion Eg(t, k) = (0(i, k) 0(i, — k)) using the r approach 
that is valid for large Peclet and Reynolds numbers. Us- 
ing Eq. (1411) written in a Fourier space we derive equation 
for the instantaneous two-point second-order correlation 
functions: 



dt 
dE e 

dt 



= <u < (k)7(-k)>+7W^ (m) (k), 
= 2{9{k)I{-k))+ME < e III \k) , 



(42) 
(43) 



where M.F^ 
72 



(HI) 



(k) 



= ~[{<hQ) + ((dui/dt)6) - 
D{uiV*9)] k and ME [ e UI \k) = -[(9Q) - D(9V 2 9)] k 
are the third-order moment terms appearing due to the 
nonlinear terms which include also molecular diffusion 
term. 

The equation for the second moment includes the first- 
order spatial differential operators Ai applied to the 
third-order moments F^ 111 ^. A problem arises how to 
close the system, i.e., how to express the third-order 
terms A4F^ IIT ' through the lower moments F^ 11 ' (see, 
e.g., [l^-dll). We use the spectral t approximation which 
postulates that the deviations of the third- moment terms, 
M.F i < III }(k), from the contributions to these terms af- 
forded by the background turbulence, AiF^ 111 ' ' (k), can 
be expressed through the similar deviations of the second 
moments, F^k) - F [IIfl \k): 

MF iHI \k)-MF (III <°){k) 
1 



r r (k) 



F^ H) {k)~ F^ n ^{k) , (44) 



(see, e.g., [88H90J]), where r r (k) is the scale-dependent 
relaxation time, which can be identified with the corre- 
lation time r(k) of the turbulent velocity field for large 
Reynolds and Peclet numbers. The functions with the 
superscript (0) correspond to the background turbulence 
with a zero gradient of the mean temperature. Validation 
of the r approximation for different situations has been 
performed in numerous numerical simulations and ana- 
lytical studies (see, e.g., review [9l|; and also discussion 
in [H, Sec. 6). 

Note that the contributions of the terms with the su- 
perscript (0) vanish because when the gradient of the 
mean temperature is zero, the turbulent heat flux and the 
temperature fluctuations vanish. Consequently, Eq. (|44[) 

for MF < > III) {k) reduces to MF i l III \k) = ~F l {k)/ T (k) 
and ME g IH \k) = -£^(k)/r(fc). We also assume that 
the characteristic time of variation of the second mo- 
ments Fi (k) and Eq (k) are substantially larger than the 
correlation time r(k) for all turbulence scales. Therefore, 
in a steady-state Eqs. (l4"2"j) and (1431) yield the following 



formulae for the functions Eg(k) and Fj(k): 

<0(k)0(-k)) = 2r 2 (fc)(^(k)%(-k))(V l T)(V J T), 

(45) 

(t*(k)0(-k)) = -TWiuiMu^-kyVjT . (46) 

Now we use the Kolmogorov model for the turbulent cor- 
relation time, r(fc) = 2ro (fc/fco) -2 ^ 3 , where tq = £q/uq 
is the characteristic turbulent time, the wave number 
k = 1/£q, the length £ is the maximum scale of ran- 
dom motions and uq is the characteristic velocity in the 
maximum scale of random motions. We also take into 
account that (ui(k)0(-k)) = 7~ 1 (u i (k) 0(-k)). Substi- 
tuting Eqs. (gS]) and d46j) into Eqs. (01]) and ((MJ), respec- 
tively, and using Eqs. (|2"T1) and (|35p we arrive at the fol- 
lowing expressions for the functions B(R) and U (S) (R) : 



fl(R) = -6/»t Aw zz (R), 



(47) 
(48) 



where ii zz (R) = (w z (x)it z (x + R)). In derivation of 
Eqs. (|47|) and (|48|) we take into account the following 
model for the second moments of turbulent velocity field, 
(ui(k) Uj(—k)) in k space: 



i(k) Uj (-k)) = 



<E{k) 
8nk 2 



k 2 



(49) 



where the energy spectrum function is E(k) — 
(2/3) k^ 1 (k/ko)~ 5/3 . In derivation of Eq. g7]) we also 
used the following identity: 



J r 3 (k) k 4 E(k) exp (ik-R) dk = -8t A E(R) 



(50) 



The theoretical study performed in this subsection allows 
us to determine the source function I(R) determined by 
Eq. CCT 



D. Equation for two-point correlation function of 
particle number density 

The theoretical analysis performed in previous subsec- 
tion has shown that the source function I(R) determined 
by Eq. (|12p . can be rewritten in the following form: 



J(R) 



Dzz(R) In 2 Re-6r 2 Au zz (R) 



N 2 

TO ' 



4a 2 



3 In Re 



(51) 
(52) 



where u zz (R) = (u z (x) w z (x + R)) and we have consid- 
ered the case a 2 ^> 1. In derivation of Eq. (fSTj) we used 
Eqs. (gUl and (gg ]) . 

Let us discuss the assumptions underlying the em- 
ployed model of particle transport in turbulent flow. We 
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use the tensor of turbulent diffusion (R) for isotropic 
and homogeneous turbulent flow. In our experiments the 
velocity field is weakly anisotropic. A main contribution 
to the magnitude of fluctuations of particle number den- 
sity is due to the mode with the minimum damping rate 
[63 . |85| . This mode is an isotropic solution of Eq. ([6]). 
Consequently, it is plausible to neglect the anisotropic 
effects. This assumption is also supported by our mea- 
surements of the two-point second-order correlation func- 
tion $(t,R). The turbulence parameters £q and uq vary 
slowly in the probed region. 

Note that the mechanism of mixing related to the tan- 
gling of the gradient of the mean particle number den- 
sity is quite robust. The properties of the tangling are 
not very sensitive to the exponent of the energy spec- 
trum of the background turbulence. The requirements 
that turbulence should be isotropic, homogeneous, and 
should have a very long inertial range (a fully developed 
turbulence), are not necessary for the tangling mecha- 
nism. Anisotropy effects can complicate the theoretical 
analysis, but do not introduce new physics in the cluster- 
ing process. The reason is that the main contribution to 
the tangling clustering is at Kolmogorov (viscous) scale 
of turbulent motions. At this scale turbulence can be 
considered as nearly isotropic, while anisotropy effects 
can be essential in the vicinity of the maximum scale of 
turbulent motions. 

Using these arguments, we consider the tensor D^R) 
in the following form: 



D„(R) = D T [F(R) + F c (R)]S i: 



, RjRj 

~~r 2 ~ 



RF' ( KiKj \ 



RiRi 



(53) 



where D T = £ Q w /3, F(0) = 1 - F c (0) and F' = dF/dR. 
The function F C (R) describes the compressible (poten- 
tial) component, whereas F(R) corresponds to vortical 
(incompressible) part of the turbulent diffusion tensor. 

Now let us study a zero-mode (i.e., a mode with 
d$/dt = 0). Using Eqs. ([6]) and (|53l) we derive equa- 
tion for the two-point second-order correlation function 
$(i?) written in a dimensionless form: 



1 



M{R) 



-/(i?){54) 



where time t is measured in units of To = £o/uq, distance 
R is measured in units of £q, and 



1 



M(R) 
x(R) 
I(R) 



2 2 
Pe + 3 
M(R) 



[1 - F - (RF C )'] , 
(F - 2F C )' , 



In 2 Re (l - 
2 A (F U (R) + 



3 Rx(R) \ 

2M(R) M(R) ) 
R K( R )' 



(55) 
(56) 

(57) 



and Pe = uo£o/D m is the Peclet number and F U (R) is 
the longitudinal correlation function of particle velocity 
field. In derivation of Eq. (f54| we have taken into account 
that for large a the term \U (A \R)^'\ < \B(R) $| for all 
scales. 

The two-point correlation function $(R) satisfies the 
following boundary conditions: $'(i? = 0) = and 
$(i? — ¥ oo) = 0. This function has a global maximum at 
R = and therefore it satisfies the conditions: 

$"(i? = 0)<0, $(i? = 0) > \<f>(R > 0)| . 

Particular formulas for the functions M(R), x(R) an d 
I(R) depend on the functions F(R), F C (R) and F U (R). 
For instance, we may choose these functions in the fol- 
lowing form: 



F(R) 
Fc(R) 

f(R) 
F U (R) 



1 



1 +<7 7 

l + a n 



cxp[-/(i?)] , 
exp[-a c /(i?)] 



R 2 



Re-i/a + R 2/3 ' 
R 2 



exp 



Re- 1/2 + R^ 3 



(58) 
(59) 

(60) 

(61) 



Equation (|60[) is similar to the interpolation formula de- 
rived by Batchelor for the correlation function of the 
velocity field that is valid for a turbulence with Kol- 
mogorov spectrum in the inertial range and for random 
motions in the viscous range of scales (see, e.g., [86l. [87| ) . 
In particular, in the inertial range of turbulent scales, 
Re~ 3 / 4 <C R <C 1, the correlation function for the tur- 
bulent diffusion tensor is F(R) cx 1 — i? 4 / 3 , where R is 
measured in units of maximum scale of turbulent mo- 
tions £q. The corresponding correlation function for the 
turbulent velocity field F U (R) cx 1 — i? 2 / 3 . The difference 
between the scalings of the turbulent diffusion tensor and 
the correlation function of the turbulent velocity field is 
caused by the scaling of correlation time r(i?) cx R 2 / 3 . 
On the other hand, in the viscous range, R <§; Re -3 ^ 4 , 
the correlation function for the turbulent diffusion tensor 
is F(R) cx 1 — Re 1 / 2 /? 2 . This correlation function is sim- 
ilar to that for the velocity field because in the viscous 
range the correlation time is independent of scale. On 
the other hand, for large scales R ^ 1 there is no turbu- 
lence, so that for R > 1 the functions F(R) and F C (R) 
should sharply decrease to zero. 

The particular choice of the functions F(R), F C (R), 
F U (R) and f(R) in the paper describes the well-known 
properties of turbulent velocity field obtained from lab- 
oratory experiments, numerical simulations and theoret- 
ical studies (see, e.g., [H, H3])- In this study we choose 
the exponential form for the functions F(R), F C (R) and 
F U (R). The final results are not sensitive to the form 
for these functions at large scales R ^> 1. Note also that 
the particular choice of the function f(R) is also not very 
important. It should describe turbulent motions in the 
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inertial range (e.g., a turbulence with Kolmogorov spec- 
trum in the inertial range) and random motions in the 
viscous range of scales. The parameter a c in Eq. (|59|) 
characterizes the different spatial scalings of the com- 
pressible and incompressible parts of the tensor of the 
scale-dependent turbulent diffusion. Our analysis has 
shown that the two-point second-order correlation func- 
tion of the particle number density is weakly dependent 
on the parameter a c . 

Analysis performed by Kraichnan in 80] for the delta- 
correlated in time turbulent velocity field, showed that 
the characteristic damping rate of the particle number 
density fluctuations is very high 7 ~ t^ 1 . The latter 
implies that the level of these fluctuations is very low 
~ (7 To) 1 ~ Re -1 / 2 . In a real flow with a finite corre- 
lation time the degree of compressibility of the turbulent 
diffusion tensor a T ^ 0, the characteristic damping rate 
of the fluctuations of the particle number density is not 
high 7 < Tq 1 , and, therefore, the level of these fluctua- 
tions is not small. 



E. Asymptotic analysis and numerical solution 

Let us perform the detailed asymptotic analysis of 
the solution of Eq. (|54|) for the two-point second-order 
correlation function $(i?) in different ranges of scales. 
There are several characteristic regions in the solution 
for the correlation function $(i?): (i) the viscous range 
< R < Re" 3/4 ; (ii) the inertial range Re" 3/4 < R < 1 
and (iii) the large scales R 1 whereby there is no 
turbulence. Here R is measured in units of the maxi- 
mum scale of turbulent motions £q. The dimensionless 
two-point second-order correlation function $(i?) is de- 
termined by Eq. (|54|) . Equation (|54|) in the viscous range 
a p < R < Re -3 / 4 reads: 



R 2 t>" + 2R& + B <f> = -B a 



(62) 



where we take into account that in the viscous range of 
scales the functions M(R), x(i?) and I(R) are given by 
the following formulas: 



Mr) = l^^^ 2 ) 

X(R) = - 2 -^M(R)Re 1 / 4 R, 



I(R) 



B a = 20J*Re , 

2(7-1 
1 + <T„ 



1 + 3(7 



3(1 



(63) 

(64) 
(65) 
(66) 



and a = a c a T . Here we use asymptotics of the functions 
F(R), F C {R) and F U (R) determined for R < Re" 3/4 [see 
Eqs. (|55 ]I -([dT 1) ]. The solution of Eq. reads 



•W)-A(£) 



-1/2 



R 



1 . (67) 



The condition $'(i? = 0) = yields ipi = 
— arctan(l/2y / So) + irk. 

In the inertial range Rc~ 3/4 < R < 1, Eq. dM} reads 

y 2 &'(y) + 4y <l'(y) + /3%) = , (68) 

where y = i? 1 / 3 , the function $ = <f> + 1, (3 = 
220/*/(9/3m) and we take into account that in the iner- 
tial range the functions M(R), x{R) an d FyR) are given 
by the following formulas: 



1 



M{R) 

x(R) 

I(R) 



= PmR a/ \ 

= _Mx M (i?)i? 4 / 3 , 

9 

81 

2(3 + 7cr) 



9(l + cr T ) ■ 
Equation (|68)l has the following solution: 



(69) 
(70) 
(71) 
(72) 



/R\~ 1 / 3 r 1- 
$(i?) = A 2 (-J cos[v^ 



i?\-V3 



V>2 



(73) 



where £ v is the Kolmogorov length scale. For large scales, 
R > 1, there is no turbulence. The condition $(R) —> 
at R — > 1 yields: A2 ~ Re 1 / 4 . Matching the functions 
<I>(i?) and $'(i?) at the boundary of the above-mentioned 
regions, i.e., at R ~ £ n yields ip2 ~ tt/2 — \fP + 7rm. 

Now we solve Eq. (|54j) numerically in order to de- 
termine the two-point second-order correlation function 
$(i?). The normalized correlation functions of particle 
number density for different parameters a T are shown in 
Figs. [2] and |U The numerical solution for the correlation 
function $(-R) is in agreement with the obtained exper- 
imental results and with the results of the asymptotic 
analysis. 

We would like to stress that in the above analysis we 
have neglected the corresponding contributions which de- 
scribe the pure inertial clustering (i.e., clustering without 
imposed mean temperature gradient), because we study 
the case a 2 3> 1, i.e., we consider the case of a small 
Stokes number (that in our experiments is St = r s /r^ = 
5.9 x 10 -2 ) and not very large Reynolds numbers (that 
in our experiments is Re = uq£qIv — 250). In this case 
contributions which determine the tangling clustering are 
much larger than that of the pure inertial clustering. 

Indeed, the main contribution to the inertial clustering 
determined by the function B(R) in the scales a p < R < 
£ v , can be estimated as 



B, 



inertial 



7 ^-Re 1 / 2 ^20a u Re 1 / 2 

(1 + Ou) 

2(4 Re 3 / 2 , 
c 4 In 2 Re 



(74) 
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FIG. 4: Normalized two-point second-order correlation func- 
tion determined from our model for a c = 1 and different 
values of the parameter a T : a T = 0.5 (solid line), a T — 1 
(dashed line), a T — 1.5 (dashed-dotted line) and a T — 2 
(dotted line). Normalized two-point second-order correlation 
function $(i?) determined in our experiments (filled squares 
as in Fig. 0). 



[see Eq. (42) in [42:]], where c s is the sound speed, and we 
have taken into account that a u ~ St 2 -C 1 . On the other 
hand, the main contribution to the tangling clustering in 
the scales a v < R < £ n is estimated as 



B 



tangling 



30 



rp2 



In 2 Re 



Re , 



(75) 



where we used Eqs. (|52"|) and (j6"5")) . Therefore, the ratio 
/Bt.an<riin<r is of the order of 



inertia] /^tangling 
-^inertial 



B 



Re 1 / 2 ^ 



tangling 



4 (V,T) 



(76) 



For the parameters pertinent to our experiments, the ra- 
tio /3* ~ 5 x 10 -9 is negligibly small. This is the reason 
why we have neglected the inertial clustering effect in 
our study. Note that for atmospheric turbulence with a 
strong temperature inversion (with a mean temperature 
gradient of the order of 1 K per 100 m) the parame- 
ter /3* ~ 10 -2 , while for atmospheric turbulence with 
a mean temperature gradient of the order of 1 K per 
1000 m, the parameter /3* ~ 1. In these estimates we 
take into account that the characteristic parameters of 
the atmospheric turbulent boundary layer are: integral 
(maximum) scale of turbulent flow £q ~ 10 — 100 m; the 
turbulent velocity in the integral scale uq ~ 0.3 — 1 m/s; 
Reynolds number Re - 10 6 - 10 7 (see, e.g., (92|.[93T|). 

The small yet finite clustering effect (see unfilled circles 
in Fig. [3]) which has been observed in our experiments in 
the absence of the mean temperature gradient can be ex- 
plained by sedimentation of particles in the gravity field. 
In this case the steady state solution of the equation for 
the mean particle number density reads: 



W 9 iV- 



L> ViV = 



(77) 



where W s is the terminal fall velocity of particles. This 
equation yields the following mean particle number den- 
sity profile for a homogeneous turbulence: 



N{z) = N exp 



z 



(78) 



with the characteristic scale of the mean particle num- 
ber density, Ljy = D T /W g ~ 10 cm for 10 fim parti- 
cles. Tangling of a gradient of the mean particle number 
density by velocity fluctuations causes particle clustering 
due to a combined effect of gravitational-tangling cluster- 
ing. However, in the turbulence with the imposed mean 
temperature gradient the effect of tangling clustering is 
considerably stronger (see filled squares in Fig. [3]). 

Indeed, the dimensionless equation for the two- 
point second-order correlation function ^(i?) for the 
gravitational-tangling clustering reads: 



1 



M(R) 



$" + 2 (i + xCR))<I>'] =-I g (R), (79) 



where the functions M(R) and x(R) are determined by 
Eqs. J55J and and 



1 - 



Rx(R) 



r 2 

' 



2M(R) M(R) 
9VF 2 



4L% 4«g 



(80) 
(81) 



The solution of Eq. (f79")) in different ranges of scales is as 
follows: 



HR) 



To 



R R 



2/?m Re 1/2 

for the viscous range of scales a p < R < Re~ 3 ^ 4 
97 



(82) 



&{R) - B 1 R- K 



2(3k + 2) 



fl 2 ' 3 , K =-L_^,(83) 



3 + 7cr 



1: 



for the inertial range of scales Re" 3/4 < R ■ 
and $(i?) = -B 2 /R for R > 1, where B 1 ~ 
3-To lnRe/4Re (3K+2)/4 and B 2 ~ 9/ /2(3k + 2). The cor- 
relation function for the gravitational-tangling clustering 
at very small scales is $(a„) ~ (Jo/v^Re) \n(£o/a p ), that 
in our experiments is of the order of ~ 0.1. The latter 
value is much smaller than that for the tangling cluster- 
ing in the turbulence with the imposed mean temperature 
gradient. 



V. DISCUSSION AND CONCLUSIONS 

In the present experimental and theoretical study we 
have found a new type of particle clustering, tangling 
clustering of inertial particles, that occurs in a stably 
stratified turbulence with imposed mean vertical tem- 
perature gradient. Fluctuations of the particle number 
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density are generated by tangling of the large-scale gra- 
dient, VAT, by velocity fluctuations. The gradient of the 
mean particle number density is formed in the stratified 
turbulence with imposed mean vertical temperature gra- 
dient due to the phenomenon of turbulent thermal diffu- 
sion. The tangling clustering of inertial particles is also 
enhanced by the additional tangling of the mean tem- 
perature gradient by velocity fluctuations, that results in 
generation of the temperature fluctuations. Therefore, 
the heat flux due to the imposed mean temperature gra- 
dient, plays a twofold role, i.e., it causes formation of 

(i) the gradient of the mean particle number density and 

(ii) a non-zero two-point correlation function of the diver- 
gence of the particle velocity field, B(R) = (t&(x) fo(y)), 
where b = div v. The latter is the main source of the 
particle clustering in small scales. 

There are two contributions to the correlation function 
-B(R). The first contribution is caused by particle iner- 
tia, so that the particles inside the turbulent eddies are 
carried out to the boundary regions between the eddies 
by the inertial forces. This is the main mechanism of the 
inertial clustering in isothermal turbulence. The second 
contribution to the correlation function -B(R) is due to 
the temperature fluctuations produced by tangling of the 
mean temperature gradient by velocity fluctuations. The 
latter results in the tangling clustering of inertial parti- 
cles. It must be emphasized that the heat flux and parti- 
cle inertia are two important ingredients which cause the 
tangling clustering. 

We have shown that in the laboratory stratified tur- 
bulence the tangling clustering is much more stronger 
than both, the inertial clustering and the gravitational- 
tangling clustering occurring in isothermal turbulence. 
In particular, in our experiments the correlation func- 
tion for the particle clustering in isothermal turbulence 
is much smaller than that for the tangling clustering in 



non-isothermal turbulence. In the stably stratified turbu- 
lence with imposed mean temperature gradient, we have 
found particle clusters with the size that is of the order 
of several Kolmogorov length scales. 

The clustering described in our study is found for in- 
ertial particles with small Stokes numbers and with ma- 
terial density that is much larger than the fluid density. 
The Reynolds numbers, based on turbulent length scale 
and rms velocity, Re = 250, in our experiments is smaller 
than that in experiments in isothermal turbulence de- 
scribed in [H, HI, [38[ . Probably these arc the reasons 
for weak inertial clustering which has been observed in 
our experiments. It must be emphasized that inertial 
and tangling clusterings in our study are understood in 
the statistical sense, i.e., we employ an ensemble aver- 
aging over many images with the instantaneous particle 
distributions rather than analyze a single instantaneous 
image. 

In the present study we have developed a theory of 
the tangling clustering of inertial particles that is based 
on the analysis of the two-point second-order correlation 
function of the particle number density. The theoretical 
predictions are in a good agreement with the obtained 
experimental results. 
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